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ABSTRACT 


The  fundamental  mechanisms  by  which  erosion  occurs  in  brittle  solids  have  been 
examined. 

The  erosion  process  in  brittle  solids  is  studied  by  considering  the  extent  to  which  a 
single  impacting  particle  produces  fracture  in  the  surface  of  a  brittle  solid.  By 
drawing  on  recent  work  on  the  location  of  fracture  in  brittle  solids,  it  is  possible 
to  predict  the  effect  of  such  variables  as  eroding  particle  size  and  velocity,  angle 
of  impact,  material  properties,  and  residual  stress.  Also,  in  the  case  of  glass 
it  is  possible  to  predict,  approximately,  the  particle  size  at  which  a  transition 
occurs  from  brittle  to  ductile  behavior. 
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INTRODUCTION 


The  erosion  of  a  surface  by  a  stream  of  solid  particles  has  received  considerable 
attention  in  the  past  decade.  About  ten  years  ago,  the  primary  motivation  for 
erosion  studies  was  the  severe  erosive  wear  that  occurred  in  the  equipment  used 
for  the  catalytic  clacking  of  oil.^)  Subsequently  other  situations  have  arisen  in 
which  erosion  has  been  a  problem;  for  example,  in  rocket  nozzles  and  most 
recently  in  the  compressors  of  helicopter  engines.  The  economic  importance  of 
erosion  in  these  and  other  applications  had  led  to  many  valuable  experimental 
studies.  However,  a  detailed  understanding  of  any  type  of  wear,  such  that  predic¬ 
tions  can  be  made,  depends  on  a  knowledge  of  the  mechanism  by  which  material  i® 
being  removed.  Unfortunately,  in  many  erosioq'  studies  this  aspect  has  been  treat¬ 
ed  superficially,  if  at  all.  1 

Clearly,  the  mechanism  of  erosion  will  noc  be  the  same  for  different  types  of 
materials,  and  to  illustrate  this,  Figure  1  shows  typical  data  for  ductile  and 
brittle  materials.  A  ductile  metal  can  undergo  extensive  plastic  deformation,  so 
during  erosion  by  rigid  abrasive  grains  we  might  expect  the  removal  of  material 
by  chip  formation  or  ploughing  as  in  metal  cutting  or  grinding.  By  contrast,  in  a 
brittle  solid  only  elastic  deformation  occurs  prior  to  fracture,  and  material  must 
be  removed  by  the  propagation  of  cracks  ahead  of  and  around  the  impacting  particle. 
The  concepts  of  ideally  ductile  and  ideally  brittle  behavior  are,  of  course,  over¬ 
simplified,  but  they  do  describe  to  a  close  approximation  the  behavior  of  many  real 
materials  and  allow  analytical  solutions  to  be  developed.  Certainly,  the  mechan¬ 
isms  by  which  brittle  and  ductile  solids  erode  must  be  understood  before  it  will  be 
possible  to  make  an  analysis  for  intermediate  behavior. 

In  this  report  we  will  first  examine  the  mechanism  by  which  brittle  solids  erode 
and  will  then  turn  to  the  question  of  ductile  solids.  We  consider,  initially,  only  the 
case  of  perpendicular  impact  (a  =  90  deg  in  Fig.  1)  and  later  will  extend  the  results 
to  the  more  general  case  of  oblique  impact. 


RELATIVE  VOLUME  REMOVAL 


Figure  1.  Volume  Removed  From  a  Surface  by  Erosion  aB  a  Function  of  the  Angle 
of  Impingement  of  the  Eroding  Particles  (Shown  Schematically  For 
Ideally  Ductile  and  Ideally  Brittle  Materials) . 
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EROSION  OF  BRITTLE  SOLIDS 


In  contemplating  an  analysis  of  brittle  erosion,  one  might  think  first  of  an  energy 
approach  in  which  the  volume  removal  is  taken  as  proportional  to  the  kinetic 
energy  of  the  impacting  particles  or  to  the  kinetic  energy  lost  by  the  impacting 
particles  during  erosion.  This  type  of  approach  is  often  followed  in  rock¬ 
drilling  studies  and  enables  crude  predictions  of  volume  removal  to  be  made. 
However,  the  values  reported  for  "energy  required  per  unit  volume  crushed",  i.e. 
the  specific  energy,  vary  greatly  for  a  given  type  of  rock,  and  it  seems  most 
-  unlikely  that  we  could  make  detailed  predictions  about  the  role  of  particle  size, 

f  particle  velocity,  impact  angle,  etc. ,  on  this  basis.  Studies  of  the  crushing  of 

<  brittle  solids  show  that  the  essential  energy  requirement,  that  to  produce  fresh 

surface,  is  a  very  small  and  unpredictable  fraction  of  the  total  energy  expended  in 
5  crushing1 2).  So  unless  fragment  size  and  surface  area  can  be  predicted  in 

advance,  it  would  appear  to  be  most  unrealistic  to  use  a  model  based  on  energy  to 
study  erosion. 

In  the  fields  of  rock  mechanics  and  coal  mining,  analyses  have  been  made  of  the 
chipping  that  occurs  under  an  indenting  tooth  or  wedge.  The  solutions  apply  only 
when  the  included  angle  of  the  wedge  is  small,  so  that  chipping  rather  than 
crushing  occurs,  and  do  not  appear  to  be  applicable  to  the  problem  we  are  con¬ 
sidering. 

Since  volume  removal  in  brittle  solids  occurs  by  fracture,  we  decided  to  approach 
erosion  by  looking  at  the  cracking  that  occurs  when  an  idealized  abrasive  particle 
of  spherical  shape  strikes  the  surface.  This  led  us  into  a  study  of  the  fracture  of 
brittle  solids  by  spherical  indenters.  Some  of  this  work  will  be  summarized  and 
applied  to  make  erosion  predictions. 


V,e  consider  the  case  of  a  sphere  pressed  against  a  semi-infinite  solid  by  a  static 
force  P  or  striking  the  plate  with  velocity  U.  The  elastic  solution  due  to  Hertz 
and  Huber  is  well  known  (e.g. ,  it  has  been  summarized  by  Timoshenko^)).  It 
shows  that  the  surface  stresses,  prior  to  fracture,  are  compressive  within  the 
contact  area  and  that  outside,  on  the  free  surface,  the  radial  stress  is  tensile 
while  the  circumferential  stress  is  of  equal  magnitude  but  compressive.  Below 
the  surface,  in  the  z-direction,  the  tensile  stress  decreases  very  rapidly. 

Figure  2  shows  the  variation  of  the  radial  tensile  stress,  on  the  free  surface, 
with  radius  R.  It  is  known  ^  that  the  impact  problem  is  pseudo  static  up  to  very 
high  velocities;  hence  The  maximum  load  P  corresponding  to  a  given  impact  velocity 
U  may  be  shown  to  be 
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where 


is  the  mass  density  of  the  sphere 
v  ^  is  the  Poisson's  ratio  of  the  sphere 
v 2  is  the  Poisson's  ratio  of  the  surface 
E^  Is  the  Young's  modulus  of  the  sphere 
Eg  is  the  Young's  modulus  of  the  surface. 


Fracture,  when  it  first  occurs,  appears  as  a  ring  on  the  surface,  usually  well 
outside  the  contact  area,  and  spreads  out  below  the  surface  to  form  the  frustum 
of  a  cone.  If  loading  is  continued,  the  initial  crack  spreads  deeper  into  the  solid 
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and  additional  concentric  ring  cracks  form  at  successively  greater  radii.*  Even¬ 
tually,  at  a  high  enough  load  or  velocity,  extensive  crushinj^occurs  under  the 
indenter.  This  behavior  is  illustrated  in  Figure  3,  where  We  have  used  a  relatively 
large  spherical  particle  for  experimental  convenience.  For  smaller  spheres,  the 
velocities  at  which  initial  cracking  anu  subsequent  crushing  occur  will  be  much 
greater  than  the  values  shown  in  Figure  3  and  also  much  greater  than  the  velocities 
encountered  in  most  applications  involving  erosion.  For  this  reason,  we  base  our 
analysis  of  erosion  on  the  type  of  cracking  shown  in  Figure  2. 

In  attempting  to  predict  ring  cracking,  there  is  little  value  in  the  classical  concept 
that  a  brittle  solid  fractures  when  the  tensile  stress  reaches  a  certain  magnitude. 
The  load  required  to  produce  ring  cracking  in  glass  shows  a  great  deal  of  scatter. 

In  addition,  there  is  a  strong  size  effect  such  that  if  the  indenter  size  is  decreased 
from,  say,  1  inch  to  0.01  inci,  there  is  a  threefold  increase  in  the  value  of  oa, 
the  stress  at  the  rim  of  the  contact  area,  computed  from  the  mean  fracture  load. 
These  observations  are  not  new,  about  30  years  ago  they  led  Weibull^)  to 
propound  his  statistical  theory  of  brittle  strength. 

Basically,  We  (bull's  idea  was  that  the  flaws  are  distributed  at  random  in  brittle 
solids.  These  flaws  are,  in  general,  too  small  to  be  detectable,  and  so  one  cannot 
relate  flaw  size  to  strength  as  in  the  familiar  "Linear  Elastic  Fracture  Mechanics" 
approach.  Rather,  one  has  to  test  a  large  number  of  specimens  and  determine  the 
probability  of  failure  for  a  given  stress  level  and  specimen  size.  As  the  size  of 
the  stressed  region  decreases,  the  probability  of  finding  a  severe  flaw  decreases 
and  as  a  consequence  the  mean  stress  at  failure  increases.  Strangely  enough, 
although  the  ring  crack  led  Weibuli  to  his  statistical  treatment  of  brittle  strength, 
he  did  not  analyze  this  problem  in  any  detail.  Recently,  the  ring -cracking  of  glass 
has  been  studied  in  some  detail^),  and  it  has  been  shown  that  the  mean  fracture 
load,  the  standard  deviation  of  fracture  load,  and  the  mean  location  of  the  initial 
ring-crack  for  indenters  of  various  radii  could  be  predicted  from  bending  test 
data.  Here,  to  illustrate  the  procedures  followed,  and  to  obtain  a  result  which  we 
will  use  later,  we  compute  the  distribution  of  the  stress  oa  at  the  rim  of  the  con¬ 
tact  area. 

Following  Weibuli^,  we  take  the  probability  F(a)  that  a  specimen  of  volume  V 
fractures  under  a  uniform  tensile  stress  <oas 


*  We  have  observed  this  type  of  cracking  during  tests  on  glass  and  by  sectioning 
aluminum  oxide  specimens  after  testing. 
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Figure  3.  Variation  of  Approach  of  Center  of  1/8  Inch  Steel  Ball  to  Undeformed 
Region  of  Glass  Plate  With  Load.  The  Impact  Velocities  Correspond¬ 
ing  to  the  Applied  Loads  Were  Calculated  From  Equation  1. 


a<au 


a  >au 


where  ou,  aQ,  m  —  the  parameters  of  the  probability  distribution  —  are  assumed 
to  be  constants  for  a  given  material.  If  these  constants  are  estimated,  e.g. ,  a 
number  of  tension  tests  on  specimens  of  a  given  size,  then  this  equation  may  be 
used  to  predict  the  probability  of  failure  at  a  given  level  for  other  specimen  sizes. 

When  the  tensile  stress,  a,  is  not  uniform  but  varies  over  the  volume,  then  the 
expression  (a  -  ffu/^0)m  is  integrated  over  the  volume  stressed  in  tension. 
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In  this  form  the  Weibuil  distribution  may  be  used,  for  example,  to  m^ke  predictions 
about  the  distribution  of  bending  strength  values. 

Glass  is  a  brittle  solid  which  is  convenient  for  experimental  work.  However,  un¬ 
like  polycrystalline  ceramics  in  which  the  strength-imparing  flaws  are  distributed 
throughout  the  volume,  glass  inevitably  fails  due  to  surface  flaws.  Thus,  if  the 
preceding  equations  are  replied  to  glass,  V  must  be  taken  as  the  surface  area 
stressed  in  tension. 


For  ball  indentation  on  glass,  the  Weibuil  distribution  takes  the  form 


F(„a> 


,  4 

=  1  -  e  ' 


3,2 

K  a  /P  -  a  \m 
a  u  \ 


2irpdp 


a  >  a 
a  u 


=  0 


a  <  a 
a-  u 


when  we  have  expressed  the  radial  tensile  stress, tr,  in  terms  of  the  maximum 
tensile  stress  oa  and  integrated  over  the  surface  area  vhere  tensile  stress  exists. 
From  this  distribution  we  may  calculate,  for  example,  the  average  value  of  aa  at 
fracture. 

Average  a  =  o  =  /  <jrdF(a) 
a  a  J  a  a 

After  obtaining  the  parameters  au,  oQ,  m.  In  this  case  from  bending  tests^,  the 
computation  may  be  carried  out  and  gives  the  results  shown  in  Figure  4.  The 
agreement  with  experiment  is  reasonable  over  a  wide  range  of  indenter  sizes. 
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Figure  4.  Predicted  and  Experimentally  Determined  Values  of  (?a  as  a  Function 
of  Indenter  Size. 


It  is  also  possible  to  use  a  two-parameter  WeibXl  distribution  for  the  calculations 
by  taking  <ru  =  0*.  This  may  lead  to  less  accurate  predictions  but  often  allows 
closed-form  solutions  to  be  obtained.  In  the  case  »f  ball  indentation  on  glass 
plates,  the  two-parameter  Weibull  distribution  leads  to  the  prediction 


o  =  (constant)R 


m  +  2 


(2) 


Another  quantity  that  is  of  interest  in  erosion  analysis  is  the  location  of  the  outer¬ 
most  ring  crack.  The  subject  of  fracture  location  has  seldom  been  considered  in 
fracture  studies  and  has  never  formed  part  of  the  Weibull  theory.  However,  by 
extending  the  Weibull  approach  and  looking  at  the  probability  of  failure  of  differ¬ 
ential  elements  which  compose  the  solid,  we  may  derive  the  distribution  of 
fracture  location. 


♦The  values  of  the  parameters  o0  and  m  are  now,  i:i  general,  different  from  those 
found  for  the  three-parameter  distribution. 
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For  a  given  load  in  bail  indentation,  if  we  divide  the  region  under  tension  into 
concentric  differential  annuli,  w  can  formulate  the  distribution  of  the  outermost 
ring  crack  location  as: 


outermost  ring  crack 

=  Prob. 

the  annulus  at  location 

occurs  at  a  radius  a* 

a*  fails 

. 

x  Prob. 

all  annuli  exterior 

to  this  annulus  survive 

This  formulation  has  been  carried  out  in  detail^,  and  it  is  shown  that  aside  from 
elastic  constants  of  the  materials,  the  distribution  of  the  outermost  ring  crack 
location  G(a*)  is  also  a  function  of  Weibull  parameters  m,  ou,  aQ,  the  indenter 
radius  R,  and  the  stress  level  in  the  solid,  oa.  That  is. 


G(a*)  =  f(a  ,  a  ,  m,  a  ,  R) 
u  o  a 


(3) 


Figure  5  shows  the  calculated  result  for  G(a*)  when  a  glass  plate  (m  =  3. 2,  ou  = 
5,000,  aQ  =  3,500)  is  loaded  by  a  3/32-inch-radius  indenter  up  to  a  load  of  154 
pounds.  The  agreement  between  prediction  and  experiment  is  quite  satisf acton,', 
particularly  in  view  of  the  complex  nature  of  the  multiple  type  of  fracture  being 
considered. 

SPHERICAL  PARTICLES  IMPACTING  AT  a  =  90° 

We  are  now  in  a  position  to  analyze  the  erosion  of  brittle  solids  and  consider 
first  the  case  of  perpendicular  impact  (a  =  90°)  by  spherical  particles.  In  all 
the  work,  we  assume  that  the  eroding  particle  is  elastic  and  does  not  fracture 
during  impact.  By  examining  fracture  patterns  such  as  those  shown  in  Figure  2, 
and  noting  that  the  radius  a*  is  considerably  larger  than  a,  it  3eems  reasonable 
to  make  the  volume  damaged  by  a  single  eroding  particle  Vp  as 

V  cc  (a*2)  S 
P 

where  a*2  =  J \ a*^  d  G(a*)  is  the  average  value  of  a*,  and  S  is  the  depth  of  the 
primary  ring  crack  as  shown  in  Figure  2.  Examination  of  ring  cracks  in  glass 
indicates  that  the  depth  of  the  primary  ring  crack  is  proportional  to  the  depth  X 
to  which  the  sphere  indents  the  surface.  From  Figure  3  we  see  that  the  Hertz 
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Figure  5.  Predicted  and  Experimental  Values  for  the  Cumulative  Distribution  of 
the  Outer  Ring  Crack  When  a  Glass  Plate  is  Loaded  by  a  3/32  Inch 
Indenter  Under  154  Pounds. 

equations  provide  a  reasonable  estimate  of  the  load-deflection  relation  even  after 
cracking  occurs.  Presumably  this  is  due  to  the  fact  that  the  stresses  are  primar¬ 
ily  compressive,  and  so  even  cracked  material  will  transmit  load. 


Thus  from  the  Hertz  equations,  we  find 


S  oc  y  oc 


2/3 

P 


cc 


so  Vp<*  RU'*/5  a*^".  Since  a*  is  a  function  of  m,  cru,  oQ,  R  and  <?a,  where  oa  in 
turn  depends  on  U,  we  now  have  an  expression  for  the  volume  damaged  by  each 
particle  which  is  a  function  only  of  R,  U  and  the  Weibull  parameters  (<rQ,  au,  m). 
We  examine  this  prediction,  which  requires  numerical  computation  to  evaluate 
the  integral,  by  varying  the  Weibull  parameter  ou,  holding  the  other  four 
quantities  (aQ)  m,  R,  U)  constant.  Effectively  ctu,  the  stress  at  or  below  which 
fracture  does  not  occur,  may  be  changed  by  inducing  biaxial  compressive  stress 
in  the  surface.  This  was  done  by  shrink -fitting  steel  rings  onto  glass  disks. 

The  predicted  and  observed  values  for  relative  volume  removal  are  shown  in 
Figure  6  and  agree  very  well.  As  another  test  of  our  formulation,  we  took 
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INDUCED  COMPRESSIVE  STRESS,  x  I03 psi 


Figure  6.  Influence  of  Biaxial  Compressive  Stress  on  Erosion  of  Glass  at  a  =  90° 
by  0. 0111  Inch  Diameter  Steel  Shot  at  a  Velocity  of  270  ft/sec. 

erosion  data  reported  by  Sheldon  and  Finnie^  for  three  brittle  solids  which  were 
eroded  by  spherical  steel  shot  of  various  sizes  at  different  velocities.  The  tester 
used  in  all  our  erosion  studies  is  described  in  the  Appendix.  The  observed  volume 
5  removal  per  particle  is  plotted  against  RU^/5  (a*^)  in  Figures  7,  8,  and  9  for 

glass,  magnesia  and  graphite  respectively.  For  magnesia  and  graphite,  (a*2) 
was  computed  by  integrating  over  the  volume  stressed  in  tension  rather  than  using 
;  area  stressed  in  tension  as  in  the  case  of  glass.  If  the  volume  removed  per 

>  particle  is  linearly  proportional  to  the  predicted  volume  damaged  per  particle, 

we  expect  the  data  points  to  fall  on  a  straight  line  of  slope  unity.  For  magnesia 
[  and  graphite  this  relation  is  indeed  observed,  while  for  glass  the  points  fall  on  a 

j  straight  line  of  slope  1.3.  We  speculate  that  the  discrepancy  in  the  case  of  glass 

j  is  primarily  due  to  the  fact  that  the  strength  of  glass  in  dynamic  tests  (erosion)  is 

|  higher  than  in  static  tests  (ring -cracking  experiments).  Possibly  the  use  of  the 

Weibull  parameters  obtained  from  dynamic  tests  would  improve  the  prediction  for 
glass. 

SPHERICAL  PARTICLES  STRIKING  AT  OBLIQUE  ANGLES 

The  extension  of  the  preceding  formulations  to  the  case  of  oblique  impact  presents 
considerably  difficulty .  In  this  case  the  stress  distribution  around  the  indenter  is 
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Figure  7.  Experimentally  Observed  Values  of  Volume  Removed  Per  Particle 

(for  a  Wide  Range  of  Velocities  and  Particle  Size)  as  a  Function  of  the 
Predicted  Damage  Volume. 


Figure  8.  Experimentally  Observed  Values  of  Volume  Removed  Per  Particle 

(for  a  Wide  Range  of  Velocities  and  Particle  Size)  as  a  Function  of  the 
Predicted  Damage  Volume. 
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Figure  9.  Experimentally  Observed  Values  of  Volume  Removed  Per  Particle 

(for  a  Wide  Range  of  Velocities  and  Particle  Size)  as  a  Function  of  the 
Predicted  Damage  Volume. 

no  longer  axisymmetric  because  of  the  unidirectional  frictional  force.  Thus,  the 
ring  cracks  are  distorted,  and  at  relatively  low  values  of  horizontal  force  they 
become  "horseshoe"  shaped  (open  in  the  direction  of  sliding). 


The  stress  distribution  in  the  case  of  oblique  elastic  impact  can  be  obtained  by  a 
superposition  of  the  Hertzian  stresses  due  to  the  normal  component  of  the  indent¬ 
ing  force  on  the  stresses  due  to  the  horizontal  component.  Using  Cartesian 
coordinates  (with  xy  lying  in  the  surface  and  z  pointing  into  the  eroded  material), 
and  assuming  a  hemispherical  distribution  of  frictional  force  in  the  contact  area, 
such  that  the  total  traction  is  equal  to  the  applied  frictional  force,  Hamilton  and 
Goodman^®)  found  for  the  contribution  of  the  horizontal  component  alone  when 
sliding  takes  place  in  the  positive  x  direction: 


2  2  9 

+  v  y  )  F  +  v  (3  -  4x  r  “)H 
pop  o 


2  2-2 
2x  F  +  (1  -  4y  r  )H 
o  o 


-  2v  x“)F  +  v  (1  -  4x^  r  ^)H 
pop  o 


(4) 

(5) 

-*,6) 
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for  r  >  a  and  z  =  0,  where  v  is  Poisson's  ratio  for  the  surface  and 


1  ,  2 

2  1/2 

1  2  | 

,  2 

21/21 

=  —  a(r 

2 

-  a  ) 

+  —  r  arctan 

2  l 

a(r 

-  a  )  | 

(7) 

1  ,  2 

2  3/2 

1  4 

f  2 

2-1/2] 

1  2  2  2  1/2 

=  —  a(r 

2 

-  a  ) 

-  —  r  arctan 

4 

|a(r 

-  a  )  | 

-  —  ar  (r  -  a  ) 

(8) 

The  stresses  given  by  these  equations  may  combined  with  the  Hertzian  stresses 
°x <2xV4  ■  r'2>]  (9) 


,  i»*i*  -  r’2, 


T  = 

xy 


3  3 


s)\i 


3  -4 

a  (1  -  2i/)  x  y  r 
P 


for  r  >t,  z  =  0.  Thus,  the  stress  distribution  on  the  surface  is  completely  known 
once  the  load  P,  the  coefficient  of  friction  fi ,  and  the  contact  radius  a  are  defined. 
Stresses  for  r  <  a  and  z  <  0  are  largely  compressive  and  are  not  pertinent  in  the 
present  study  because  fracture  is  assumed  to  originate  from  a  region  in  tension. 

Since  fracture  propagation  follows  a  principal  stress  trajectory^),  it  is  assumed 
that  for  oblique  impact,  the  volume  removed  is  given  by 


V  =  (constant) 
P 


Mean  value  of  A(x  ) 
o 


where  is  the  maximum  depth  of  indentation  produced  by  the  vertical  component 
U  sin  a,  and  A(xQ)  is  the  area  bounded  by  the  outermost  principal  stress  trajectory 
(PST)  on  the  z*0  plane.  Note  that  for  zero  coefficient  of  friction  between  the 
contacting  surfaces,  the  stress  distribution  is  Hertzian,  and  the  prediction  of 
volume  removal  is  identical  with  the  formulation  for  normal  impact.  The  distri¬ 
bution  of  the  outermost  PST  can  be  derived  in  a  manner  similar  to  that  for  the 
outermost  ring  fracture;  namely,  it  is  stipulated  that  the  material  outside  the  PST 
in  question  survives,  and  the  material  within  fractures  in  at  least  one  place. 
Numerically,  however,  the  distribution  of  the  outermost  FST  is  much  more 
difficult  to  evaluate. 

Let  a  PST  be  described  by  its  intercept  on  the  negative  x-axis,  say,  x0,  as  shown 
in  Figure  10.  Thus,  in  a  region  with  multiple  cracks,  the  outermost  xD,  denoted 
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Figure  10.  Numerical  Method  Used  to  Obtain  Principal  Stress  Trajectory. 


by  Xf,  will  be  a  random  variable  whose  distribution  is  related  tc  the  properties  of 
the  two  bodies  in  contact,  and  the  amount  of  indenting  force.  Now, 

Prob.  (Xf  <  x^)  =  Prob.  (the  material  outside  the  PST  defined  by 

survives  the  given  load) 

x  Prob.  (the  material  within  fails) 


or 


-/M 


m 


dA 


Prob.  [X  <  x  1 

L  f  fj 


=  e 


all  A  2 


PST(Xf) 


=  e 


(o  ~  a  \m 

-/(M  - 

all  A  ><Wl'(x{) 


for  which  the  probability  density  function  is 
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BP 


m 


all  A  >  PST(xf) 


dA 


all  A  >  PST  (x^) 


dA 


(13) 


Obviously,  numerical  methods  are  necessary  to  evaluate  this  equation,  which 
mainly  Involves  the  integral 


all  A  >  PST(xf) 


(14) 


in  the  plane  z  =  0  and  the  region  outside  the  PST  prescribed  by  Xj  for  which  o1  > 
This  integral  involves  meticulous  calculation  and  is  a  major  task  even  with  a 
high-speed  computer.  First,  for  fi  >  0,  the  stress  distribution  is  highly  complex. 
Great  difficulty  also  arises  in  calculating  the  PST,  for  which  there  is  no  closed- 
form  solution;  then  there  is  the  evaluation  of  the  integral  itself  by  numerical 
methods.  Nevertheless,  a  program  was  prepared  to  perform  the  following 
sequence  of  calculation: 

THE  PRINCIPAL  STRESS  TRAJECTORY 


Referring  to  Figure  10,  starting  with  =  0  and  y0  =  0,  the  direction  of 
the  principal  stress  is  calculated  using  the  equations  for  the  stresses.  A  neigh¬ 
boring  point  with  coordinates  (x^',  y^  is  then  defined  by 


x  ’  =  x  +  As 
1  0 

cos  /31 

(15) 

yi*  =  yQ  +  As 

sin  Pi 

(16) 

where  As  is  an  arbitrary  increment.  At  (xj',  yj'),  the  direction  of  the  principal 
stress,  is  again  calculated,  and  a  point  supposedly  on  or  near  the  PST  is 
defined  by 


Xj'  =  x1'  +  As  sin  -/32)  sin^ 


(17) 


Y1  =  yx’  +  As  sin  (^  -  pj  cos  $2 


(18) 


j 

I 

i 


16 


This  procedure  is  then  repeated  to  obtain  another  point  (x2,  y2),  and  so  on.  The 
points  (xQ,  y0),  (xp  y^),  (x2,  y2)  .  .  .  are  then  points  on  or  near  the  PST  which 
passes  through  (Xo,y0).  It  is  obvious  that  without  Equations  (17)  and  (18),  the 
coordinates  for  the  PST,  which  would  have  been  calculated  only  by  Equations  (15) 
and  (18),  will  be  subject  to  cumulative  errors.  With  Equations  (17)  and  (18)  it 
was  possible  to  have  virtually  all  the  calculated  points  fall  on  the  PST. 

The  area  bounded  by  the  PST,  A(xc),  is  easily  found  once  the  coordinates  of  the 
points  in  the  PST  are  known.  In  fact,  A(xq)  is  found  by  summing  the  incremental 
areas  Aj,  shown  shaded  in  Figure  10,  which  were  calculated  at  the  same  time  the 
coordinates  (x{,  y.)  were  being  found.  Thus,  referring  to  Figure  10, 


Aa  = 

i 


R  +  R 
1  2 


where  8  is  in  radians,  and 


A(xo)  =  £  AA. 


ria\  ~  c\i\m 

The  integi-  l  1 1 - -)  dA  is  approximated  by  the  sum 

all  A  >1>ST(x  ) 


S2<V  ' 


X  >  X 


along 

PST(x  ) 
o 


a.  -  a 
1  u 


As  Ax 


The  summand  (Oj  -  au/aQ)m.  As  in  the  inner  sum  is  also  calculated  at  the  same 
time  that  the  coordinates  are  found.  The  implication  of  this  equation  is  that  in  the 
elemental  area  As  *  Ax,  the  principal  stress  oj  is  constant.  Equation  (21)  accounts 
for  only  the  shaded  part  of  the  total  risk  of  rupture  shown  in  Figure  11  and,  as  a 
result, will  underestimate  the  integral  by  an  amount  which  depends  upon  how  the 
PSTs  '’diverge".  At  any  rate,  the  errors  in  the  preceding  equation,  as  in 
Equations  (15),  (16),  (17),  and  (18),  can  be  minimized  by  making  As  and  Ax  as 
small  as  practical.  It  was  found  that  for  As  =  0.05x  and  Ax  =  0.01a,  the  errors 
are  insignificant  in  the  sense  that  by  making  As  or  Ax  smaller,  there  is  very 
little  change  in  the  predicted  erosion  rate. 
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Figure  11.  Error  in  the  Numerical  Method  Used  to  (  lculate  Erosion  Rate 
Under  Oblique  Impact  Conditions. 

Finally,  with  A^)  and  S  (xq)  stored  tn  the  memory  of  the  computer  for  all 
a  <  at  an  interval  of  Ax  =  0.01a,  where  XqU  is  the  value  of  beyond  which 

al  <  °u’  average  volume  removal  per  particle  is  calculated  by  taking  the  pro¬ 
duct  of  the  depth  of  indentation  and  the  average  value  of  Afx^ .  That  is. 


V 

P 


const  X 
m 


A(x  )  f  (x  )  dx 
o  Xj.  o  o 


(22) 


where  fxj(x0)  may  b®  found  once  the  value  of  S2(x0)  Is  known.  It  is  to  be  noted  that 
in  the  evaluation,  the  equation  for  S2(xQ),  a  surface  distribution  of  flaws,  is 
assumed  —  that  flaws  exist  only  on  the  surface,  or  at  least  that  those  underneath 
are  insignificant  compared  to  those  on  the  surface  insofar  as  fracture  initiation 
is  concerned.  This  assumption  is  probably  justified  in  the  case  of  glass,  where 
from  previous  experience  failure  is  known  to  initiate  invariably  from  a  surface 
flaw.  In  other  brittle  materials  such  as  alumina,  magnesia  and  graphite,  the 
location  of  fracture  is  not  clear,  and  hence  there  is  some  reservation  regarding 
this  assumption.  In  any  case,  the  principal  stress  decays  rapidly  with  the 
distance  from  the  surface  and  becomes  compressive  at  a  distance  of  about  0.05a. 
Thus,  treating  the  problem  in  this  manner  is  probably  not  too  far  from  the  actual 
condition.  The  only  unknown  parameter  in  the  analysis  is  the  coefficient  of 
friction  between  the  sphere  and  the  surface,  for  this  determines  the  horizontal 
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Figure  12.  Predicted  and  Observed  Values  of  the  Relative  Volume  Removal  as  a 
Function  of  Angle  for  Glass  Eroded  by  0. 0165 -Inch -Diameter  Steel 
Shot  at  200  ft/sec.  Predictions  are  Made  for  Three  Values  of  n,  the 
Coefficient  of  Friction.  The  Inserted  Sketch  Shows  the  Contact  Circle 
and  the  Horseshoe  Cracks. 

force.  For  impact  at  high  angles  (near  90  deg),  the  spheres  should  roll  rather  than 
slide  and  the  coefficient  of  friction  will  be  very  low.  In  Figure  12,  the  predictions 
are  compared  with  experiment,  and  it  is  seen  that  good  agreement  is  obtained  at 
angles  near  90  deg  with  a  coefficient  of  friction  of  0. 03. 

AnGULAR  ABRASIVE  PARTICLES  IMPACTING  AT  a  «  99  DEGREES 

r 

In  practice,  erosion  is  produced  by  angular  abrasive  particles  rather  than  spheres. 
Since  a  detailed  analysis  of  this  problem,  similar  to  that  developed  for  spheres, 
does  not  appear  possible  at  the  present  time,  we  present  an  approximate  solution 
following  the  approach  suggested  by  Sheldon  and  Finnie^  which  was  later  develop¬ 
ed  by  H.  L.  Oh  and  Finnie'11^. 

Surfaces  struck  by  eroding  particles  typically  have  the  appearance  shown  in 
Figure  13.  There  is  relatively  little  damage  under  the  particle,  perhaps  because 
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Figure  13.  Schematic  View  of  Cracking  Due  to  Particle  Impact . 


of  large  hydrostatic  compressive  stresses  at  this  point.  On  the  surface  surround¬ 
ing  the  indentation,  a  fairly  extensive  region  of  material  is  cracked.  Based  on 
this  observation,  we  assume  that  t}ie  volume  removed  by  an  eroding  particle, 

V  is  proportional  to  the  damaged  yplume  enclosed  in  the  spherical  cap  of  radius 
R*  and  depth  S.  That  is. 


trurgi'i  n. 

,if9ffi991J3£  bo 


2 

V;  oc  S  R* 

v 


.so. 

In  case  the  particle  is  not  a  sphere  but  is  Irregular  in  shape,  we  may  invoke 
two  radii:  radlu^^^igj}  cc^vesponda  to  a  sphere  weighing  the  same  as  the 
particle,  and  radiusTt*,  wEicFTs  effective  in  producing  fracture. 


'ri;ET  39bitl£q  “ 

Now,  we  have  to.,gg£jjg^g  ^gpd  R*;  a3  before,  we  shall  take  S  as  proportional  to 
the  depth  of  indent  trained  from  the  Hertz  solution.  We  consider  R*  to  be 

the  radius  of  a  much  larg^lr  indenter  than  the  actual  particle,  which  on  the  average 
produces  a  ring  crack  at  radius  a*.  This  assumption  is  perhaps  not  too  unreason¬ 
able  because  the  cracked  material  is  capable  of  transmitting  load.  Since  we  have 
seen  from  Figure  3  that  the  Hertz  equations  continue  to  give  a  reasonable  estimate 
of  load-deflection  behavior  even  after  considerable  cracking  has  occurred,  we 
write 
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S  «  X  «  R*  a 

a 

In  order  to  get  a  closed  solution,  we  now  use  the  two-parameter  Weibull  distribution 
(<ru  =  0)  and  from  Equation  (2)  obtain 


a  X  «  R* 
a* 


2 

m  +  2 


Eliminating^*,  we  arrive  at  a  relation  between  R*  andX: 


R*  cc  X 


m  +  2 
m  -  2 


The  depth  of  indentation  X  may  be  estimated  readily  by  assuming  that  the  Hertz 
solution  holds  true  even  after  cracking  occurs.  Then  from  Timoshenko'3 


x  «  [r6/rOiA  «4/5 


where  U  is  the  particle  velocity.  Finally,  volume  removed  per  particle  is 


3m  -  2 

V  a  S2  R*  «  X2  R*  a  X  m  "  2 
P 


6,1/5  \  m  -  2 


/3m  -  2 


[£] 


ff)  .A 


u 


4/5 


3m  -  2  \ 
m  -  2  > 


We  distinguish  two  cases: 


spherical  particle: 


R'  =  R,  V  «  R 
P 


/3m  -  2\ 
J  m  -  2  / 


U 


4/5 


/3m  -  2\ 
'  m  -  2  / 


/3m  -  2\  /3m  -  2\ 

irregular  particle:  R'  «  constant,  V  «  R6//5  m  “ 2  u4^5  m  "  2 


These  expressions  for  the  dependence  of  volume  removal  on  particle  size  and 
velocity  have  been  confirmed.  In  experiments  on  six  brittle  solids  over  a  wide 
range  of  velocities  and  particle  sizes, <i2)  it  was  indeed  found  that  the  volume 
removal  per  particle  could  be  written  as 
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for  spherical  steel  shot 


P 

a'  b 

V  oc  R‘  U  for  angular  silicon  carbide  -particles 
P 

A  comparison  of  the  predicted  and  observed  values  of  a,  a’  and  b  is  given  in 
Figure  14.  Except  for  glass  eroded  by  spherical  steel  shot,  the  correlation 
between  theory  and  experiment  is  quite  good,  particularly  in  view  of  the  number  of 
assumptions  involved. 

THE  BRITTLE -DUCTILE  TRANSITION 

An  interesting  consequence  of  the  size-effect  on  strength  of  brittle  solids  is  that 
if  the  solids  are  loaded  on  smaller  and  smaller  regions,  the  stress  required  to 
produce  fracture  may  eventually  exceed  that  required  for  yield.  Thus,  we  can 
understand  the  familiar  observation  that  microhardness  or  scratch-hardness  tests 
on  brittle  solids  may  produce  flow  rather  than  fracture. 

To  estimate  the  particle  size  at  which  ductile  behavior  should  be  noticed  in  erosion 
testing  nominally  brittle  solids,  we  need  only  the  parameters  of  the  Weibull 
distribution  and  indentation  hardness  measurements.  For  glass,  this  estimate  can 
be  made  even  more  directly  from  the  data  shown  in  Figure  4.  The  mean  value  of 
03  (the  stress  at  the  rim  of  the  contact  area)  required  for  fracture  is  related  to 
the  average  stress  p  under  the  indenter  through  the  Hertz  equations  as 

p  =  °a  2/(1  -  2v)  «  4oa  v  _ ;  (23) 

-3 

So,  for  an  indenter  of  radius  10”  inches,  we  would  estimate  oa  =  250,000  psi 
and  p  =  106  psi  from  Figure  4.  Values  of  the  flow  pressure  P(y ieid)  *or  6^a8S 
have  been  given  by  Marsh  (13)  and  are  typically  about  800,000  psi  for  ordinary 
hardness  tests  and  about  2,000,000  psi  for  hardness  tests  at  low  temperature 
where  time -dependent  effects  are  absent.  The  points  in  Figure  4  are  of  course 
for  static  tests  and  would  be  somewhat  higher  if  dynamic  loading  had  been  used. 

So  direct  comparison  is  impossible,  but  it  appears  that  yield  should  occur  before 
fracture  for  particles  in  the  range  10”4  to  10“^  inches.  Erosion  tests  on  glass 
show  just  such  behavior.  In  Figure  15  it  is  seen  that  erosion  tests  on  glass  with 
particles  5  x  10"^  inches  and  8  x  10-4  inches  in  diameter  lead  to  typically  brittle 
behavior,  while  3.5  x  iO”4-inch-diameter  particles  show  erosion  behavior  that  is 
typical  of  ductile  metals.  Other  aspects  of  the  brittle -due'  lie  transition  in  erosion 
have  been  discussed  by  Sheldon  and  Finnie(14^. 
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ANGLE  ot  ,  degrees 

Figure  15.  Weight  Removal  Per  Gram  of  Abrasive  as  a  Function  of  Particle 

Approach  Angle  a  for  Glass  Eroded  by  Silicon  Carbide  Particles  at 
500  ft/sec. 

QUANTITATIVE  PREDICTION  OF  VOLUME  REMOVAL 

In  our  treatment  of  erosion  due  to  perpendicular  Impingement  by  solid  particles, 
we  took  the  volume  removed  by  a  single  particle  as  proportional  to  the  cracked 
volume  produced  by  a  single  particle  striking  a  smooth  surface.  Thus,  at  least 
one  erosion  test  must  be  made  on  the  chosen  material  before  quantitative  pre¬ 
dictions  can  be  made.  We  took  this  approach  because  the  intersection  of  the 
fracture  ratteras  from  many  particles  will  be  a  complex  process,  and  more  work 
needs  to  be  done  on  this  topic  before  a  complete  quantitative  treatment  of  erosion 
is  possible.  However,  It  may  be  of  value  to  outline  some  preliminary  consider¬ 
ations  in  this  regard. 

Generally;  one  would  expect  that  the  extent  of  cracking  in  any  brittle  solid  would 
depend  only  on  the  stress  state  in  the  part  prior  to  fracture  and  on  the  fracture 
toughness  of  the  material.  This  latter  quantity  could  be  expressed  by  an  effective 
surface  energy  for  fracture  under  plane  strain  conditions  G|C  in.  -lb/in. 2  or  by 
the  related  value  of  the  stress  intensity  factor  KjC  lb/in.3/2,  where 
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K2  K  2 

ic  2  rc 

G  =  - —  (1  -  V  )  -  — — -  (24> 

IC  E  '  '  E 

Normally,  the  elastic  modulus  does  not  appear  in  elastic  solutions  for  the  stresses; 
but  in  the  present  case,  because  of  the  change  of  contact  area  with  load,  we  would 
expect  the  stress  state  and  hence  the  crack  depth  in  the  solid  to  depend  on  the 
modulus  as  well  as  on  the  load  (or  velocity)  and  the  radius  of  the  indenter. 

In  our  analysis  of  erosion  by  spherical  particles,  we  took  the  depth  of  cracking  as 
proportional  to  the  depth  of  indentation.  TV  - 1  is. 


Rl/3E2/3 


where  E  is  a  weighted  average  of  the  moduli  of  indenter  and  surface  given  by  the 
Hertz  equations.  A  limitation  of  this  expression  is  that  the  surface  energy  for 
fracture  is  not  included  and  hence  the  constant  of  proportionality  must  depend 
on  Kjq  (or  Gjq)  . 

By  contrast,  Barenblatt^15),  in  discussing  the  literature  on  fracture  under  a 
cylindrical  punch,  where  the  contact  area  does  not  change  with  load,  showed 
that 


,GICE> 


Very  recently,  we  have  made  further  studies  of  the  depth,  S,  of  the  initial  ring 
crack  in  glass  loaded  by  spherical  indenters.  Rather  than  finding  S  proportional 
to  the  indentation  X,  i.e. ,  S  ~  p2' 3/R1/2,  recent  work  suggests  S  ~  p5' VR1/2. 
Hence  from  dimensional  considerations,  we  speculate  that 

5/6  5/6 

S  ~  R1/3GICPl/3E1/2  ~  R1/3  KIC2/3  E1/6 


This  new  result  will  not  change  any  of  our  previous  predictions  significantly  but 
does  suggest  that  the  coefficient  of  proportionality  between  the  volume  (a*2)X 
and  the  volume  removed  per  particle  should  involve  l/(G|c*/3  E1'2). 
Representative  values  of  Gjq  are  not  easy  to  obtain  for  brittle  solids.  However, 


25 


r 


from  the  literature,  we  obtained  the  following  approximate  values  for  the  materials 
studied  in  Figures  7,  8,  and  9. 


G  in'"ib 

E  JL_ 

,gic-1/3e-1/2> 

Observed  Ratio  of 

IC  in. 

.  2 

in. 

Damage  per  Particle 
+(a*2)X 

Glass 

15  x  io-3 

10  x  106 

3 

9 

Graphite 

5  x  KT1 

1  x  106 

r* 

0 

7 

MgO 

80  x  10"3 

30  x  106 

1 

1 

A1,C>, 

20  x 10“2 

53  x  106 

0.5 

(not  tested  with  shot) 

*  Relative  to  MgO  as  unity 


Comparison  of  the  last  two  columns  shows  that  this  method  of  predicting  volume 
removal  is  far  from  precise  but  may  be  worth  pursuing  in  future  work. 
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CONCLUSIONS 


It  has  been  shown  that  by  considering  the  region  damaged  by  a  single  eroding 
particle  we  can  predict  many  aspects  of  the  erosion  of  brittle  solids.  The  effects 
of  impingement  angle,  particle  velocity,  particle  size  and  material  properties 
are  all  included  in  the  analysis.  It.  addition  it  is  possible  to  predict  the  transition 
from  brittle  to  ductile  behavior  in  the  case  of  glass,  when  the  particles  are  small 
enough. 

At  the  present  state  of  the  analysis,  at  least  one  erosion  test  must  be  made  on  the 
material  of  interest  before  erosion  rates  can  be  predicted.  Further  work  on  the 
depth  of  cracking  under  a  single  particle  and  on  the  interaction  of  cracks  produced 
by  many  particles  would  be  desirable.  With  this  information  it  may  be  possible  to 
make  completely  quantitative  erosion  predictions  based  on  fracture  test  data. 
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APPENDIX 

EROSION  TEST  EQUIPMENT 


For  completeness,  we  include  this  description  of  our  erosion  test  apparatus.  It 
is  essentially  that  given  by  Sheldon^). 

An  erosion  tester  must  be  capable  of  accelerating  a  small  but  variable  mass  flow 
rate  of  abrasive  particles  to  a  predetermined  velocity,  allow  them  to  impinge  on 
a  test  specimen,  and  then  be  exhausted  into  a  suitable  collector.  Provision  must 
be  made  in  the  tester  to  vary  the  angle  of  attack  between  the  abrasive  particle 
approach  direction  and  the  surface  of  the  test  specimen.  To  allow  rapid  testing 
of  a  large  number  of  samples,  the  tester  must  be  simple  and  easily  operated. 

Figure  16  is  a  schematic  view  of  the  apparatus  designed  to  fulfill  these  objectives. 
Figure  17  shows  photographs  of  the  same  equipment.  The  equipment  is  designed 
to  be  used  in  two  configurations.  The  First  configuration  is  sed  for  initial 
pressure-velocity  calibration  of  the  abrasive  particles,  in  which  case  no  erosion 
sample  is  in  place.  In  the  second  configuration  the  test  specimen  is  eroded  with 
the  velocity  calibration  equipment  removed. 


HIGH  PRESSURE 
AIR 


Figure  17.  Photographs  of  Erosion  Test  Equipment 


In  the  latter  mode,  the  equipment  functions  as  follows.  Referring  to  Figure  16,  a 
measured  amount  of  abrasive  grit  is  placed  in  the  pressure-tight  container  A 
through  the  screw  top.  Container  A  is  mounted  on  rubber  shock  isolation  mounts, 
and  hence  is  capable  of  a  limited  amount  of  movement.  It  is  made  to  vibrate  at  a 
desired  rate  by  means  of  a  small  pneumatic  vibrator  B.  The  vibration  of  the  con¬ 
tainer  feeds  the  particles  at  a  controlled  rate  down  the  feed  tube  and  through  the 
flexible  connection  C  into  a  mixing  chamber  at  E.  A  small  constricting  ring  in 
the  feed  tube  insures  a  uniform  flow  of  particles  into  this  mixing  chamber.  High- 
pressure  air  is  filtered  at  D  and  enters  the  mixing  chamber  through  tube  F. .  The 
abrasive  mixed  with  the  high -pressure  air  is  then  accelerated  through  a  length  of 
tubing  G  to  the  desired  velocity  and  allowed  to  Impinge  on  the  specimen  H.  The 
specimen  is  held  firmly  to  a  vacuum  plate  J,  whose  angular  position  o  is  adjustable 
from  10  to  90  degrees.  The  axis  of  rotation  lies  on  the  specimen  surface;  hence 
the  particle  acceleration  distance  remains  constant  with  changing  a.  Tne  eroding 
operation  is  carried  on  Inside  a  sealed  container  K.  This  container  is  a  box  about 
six  inches  square  and  has  provision  for  viewing  windows  or  access  doors  on  all 
four  sides.  When  erosion  tests  are  carried  out,  the  adjustable  vacuum  plate  J  is 
fitted  in  place  of  an  access  door.  This  containing  box  is  connected  indirectly  to  a 
shop  type  vacuum  cleaner  V,  which  has  sufficient  capacity  to  maintain  a  3-  to  4- 
inch  Hg  vacuum  under  all  operating  conditions.  This  vacuum  system  would, 
probably,  not  be  usjable  in  standard  form  for  high -temperature  experiments. 

After  impacting  the  specimen,  the  abrasive'  particles  along  with  all  eroded 
materials  are  drawn  out  of  the  blasting  chamber  and  deposited  by  a  conventional 
cyclone  separator  and  are  picked  up  in  the  vacuum  cleaner. 

A  typical  test  run  would  proceed  as  follows.  The  vacuum  plate  J  is  set  at  the 
desired  impact  angle  and  locked  in  position.  The  specimen  is  accurately  weighed 
and  placed  on  the  vacuum  plate,  being  held  firmly  in  position  ie  vacuum.  The 
specimen  is  normally  1/8  inch  thick  by  1  inch  wide  but  may  be  u.  long  as  3  inches 
depending  oA  the  attack  angle  a.  An  accurately  weighed  amount  of  abrasive  is 
placed  in  container  A,  and  the  container  is  sealed.  Abrasive  charges  may  range 
from  0. 05  to  100  gm.  From  previous  tests,  the  air  regulators  Rj  and  R2  are 
set  with  gages  Gj  and  G2  to  give  the  desired  abrasive  velocity  and  feed  rate 
respectively.  Then,  with  all  access  doors  in  position,  the  vacuum  cleaner,  the 
blast  air  controlled  by  switch  Sj,  and  the  vibrating  feeder  controlled  by  switch 
S2  are  switched  on,  in  that  order,  and  the  erosion  process  begins.  The  abrasive 
particles  are  normally  fed  at  about  10  grams  per  minute,  the  abrasive  feed  rate 
being  readily  observed  through  a  window  in  feeder  A.  When  all  the  particles  have 
been  exhausted,  the  test  is  complete.  Final  weighing  of  the  test  specimen  then 
determines  the  weight  of  eroded  material.  If  the  eroded  cavity  is  found  to  be  deep 
enough  to  change  the  apparent  Impact  angle  a  significantly,  the  test  is  repeated 
using  a  reduced  abrasive  charge. 


The  vibrating  feeder  works  very  well  for  particles  larger  than  about  500  mesh 
(0.0008  Inch).  In  tests  using  particles  of  500  mesh  or  smaller,  feeder  A  is  not 
used.  Instead,  the  high-pressure  air  inlet  tube  F  is  replaced  with  a  similar  tube 
having  a  small  aspirator  tube  extending  into  its  wall.  A  small  flexible  hose  is 
connected  to  the  aspirator  tube.  The  desired  charge  of  abrasive  particles  is  spread 
in  a  thin  layer  in  a  pan,  and  the  open  end  of  the  flexible  aspirator  hose  is  passed 
over  the  particles  to  entrain  them  in  the  main  air  blast.  The  abrasive  flow  rate  is 
not  as  uniform  using  the  aspirator  as  it  was  using  the  vibrating  feeder;  however, 
the  results  obtained  were  very  satisfactory.  The  aspirator  was  not  usable  at 
blast  pressure  above  about  25  psig.  This,  however,  caused  no  difficulties,  as 
these  smaller  particles  are  accelerated  to  the  desired  velocities  by  much  lower 
bla$t  pressures. 

In  all  the  erosion  tests,  the  blast  tube  G  was  12  inches  long  by  0. 182  inch  inside 
diameter.  The  tube  material  was  304  stainless  steel.  Because  the  tube  bore 
wore  rather  quickly  and  particle  velocity  was  very  sensitive  to  bore  changes,  the 
apparatus  was  constructed  to  enable  rapid  tube  changes  to  be  made.  In  this 
application,  a  tube  of  high  hardness  ceramic  material  would  have  been  more 
suitable. 

An  importact  parameter  in  erosion  analysis  is  the  abrasive  particle  velocity  U. 

This;  quantity  cannot  be  determined  with  sufficient  accuracy  using  analytical 

methods  because  of  the  significance  of  several  unknown  parameters.  For  example, 

it  is  difficult  to  account  for  the  effects  of  impact  between  the  particles  and  the 

walls  of  tube  G  or  between  the  particles  themselves.  It  was,  therefore,  decided 

to  determine  the  particle  velocities  experimentally. 

j 

1 

The  technique  used  for  measuring  particle  velocities  was  to  photograph  the 
particles  twice  using  a  double-flash  light  source  M.  Knowing  the  interval  between' 
the  light  flashes  and  by  measuring  on  the  film  the  distance  the  particle  travels 
between  exposures,  the  average  particle  velocities  are  easily  calculated. 

i 

The  light  eource  used  was  a  commercially  available  unit,  which  has  a  variable 
flash  interval  of  5  to  100  microseconds  with  a  flash  duration  of  0. 5  microseconds. 

A  Frehnel  lens  mounted  on  the  front  of  the  source  efficiently  accumulated  the 
short -rduratlon  light.  An  opaque  tube,  0.6  inch  in  diameter  and  16  inches  long 
formed  a  light-tight  path  between  the  flash  unit  and  the  viewing  box.  The  image 
of  *he  Xenon  flash  tubes  was  focused  on  the  camera  lens  by  an  adjustment  built 

into  the  unit. 

» 

\  : 

The  camera,  shown  at  N,  was  a  conventional  Speed  Graphic  with  13 -inch  bellows 
pxteusion  and  removable  lens  board.  A  100  mm  Zeiss  Luminar  lens  mounted  on 
a  3-inch  extension  tube  provided  a  capability  of  3.5X  magnification.  The  Luminar 
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lens  was  especially  suitable  for  this  application  because  of  its  extreme  sharpness 
and  its  effective  reduction  of  curvature  of  field.  The  lens  board  of  the  camera  was 
connected  to  the  view  boss  by  an  opaque  bellows  P,  variable  in  length  from  2  to 
12  inches.  This  provided  a  light-tight  path  between  the  view  box  and  the  camera, 
yet  allowed  the  camera  to  lie  focused  on  the  abrasive  stream  at  a  fixed  camera 
bellows  extension  or  magnification. 

Conventional  4x5  sheet  film  was  used  in  the  tests.  It  was  found  that  when'tj.lm 
with  rating  ASA  400  was  used,  a  lens  opening  of  f6. 3  provided  good  exposure.  For 
trial  runs,  a  Polaroid  4x5  sheet  film  adapter  was  used  in  the  Graphic  camera 
and  was  found  to  be  very  convenient.  Polaroid  3000  series  film  shot  at  fl6  was 
used  and  provided  good  resolution,  though  not  as  good  as  the  regular  sheet  film. 

In  photographing  a  stream  of  abrasive  particles,  the  primary  problem  was  in 
determining  the  two  positions  of  the  same  particle  on  the  doubly  exposed  film.  It 
1  was  often  very  easy  to  confuse  two  separate  particles  as  being  the  same  particle 
1  in  its  two  positions.  To  reduce  the  possibility  of  this  mix-up  occunlng,  it  was 
-scessary  to  reduce  the  number  of  rebounding  particles  inthe  viewing  box  K. 

This  was  effectively  accomplished  by  removing  the  specimen  holder  J  and  allowing  ». 
the  stream  of  particles  to  pass  straight  through  the  box  into  the  precipitator.  In 
addition,  by  reducing  the  feeder  pressure,  the  particle  density  was  reduced  so  only 
a  few  particles  were  visible  in  a  single  exposure. 

*  l 

The  procedure  used  to  photograph  an  abrasive  particle  was  as  follows.  The  camera 
bellows  was  set  to  give  the  desired  magnification,  and  the  camera  was  carefully 
focused  on  the  center  line  of  tube  G.  The  double -flash  light  source  was  set  to  give 
u.  particle  travel  distance  of  about  6  times  the  particle  diameter.  The  abrasive 
flow  was  started  and  the  blast  air  pressure  was  set  at  the  desired  level.  A  photo¬ 
graph  was  then  taken  by  removing  the  light  shield  from  the  film  holder  and  flash¬ 
ing  the  light  source.  Because  of  the  precautions  taken  to  eliminate  light  leakage, 
no  shutter  was  necessary  on  the  camera. 

The  particle  travel  distances  were  accurately  measured  using  a  Zeiss  measuring 
microscope.  This  instrument  permits  viewing  through  the  negative  at  magnifica¬ 
tions  up  to  80X.  The  built-in  measuring  system  permits  an  easy  measurement  to 
be  made  between  the  centers  of  gravity  of  the  particle  images  accurate  to 
i/10. 000th  of  an  inch.  Abrasive  particles  as  small  as  500  mesh  (0.0008  inch 
diameter)  could  be  viewed  with  this  instrument. 

In  all,  about  500  photographs  were  taken  and  about  3000  separate  measurements 
were  made  to  obtain  calibration  curves  of  particle  velocity  versus  air  pressure 
for  different  sizes  of  steel  shot  and  silicon  carbide  grit. 
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